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ABSTRACT 

The NASA Lewis Research Center manages for the Department of Energy, the technology and engineering 
development of large horizontal axis wind turbines. In support of this activity each wind turbine has 
a variety of information systems used to acquire, process and analyze data. In general four categories 
of data systems, each responding to a distinct information need, can be identified. The categories 
are: Control, Technology, Engineering and Performance. 

The focus of this report is on the information that can be extracted by statistical analysis of data 
obtained from the Technology and Engineering Information Systems. These systems consist of the 

following elements: (1) sensors which measure critical parameters (e.g. , wind speed and direction, 

output power, blade loads and component vibrations; (2) remote multiplexing units (RMUs) on each wind 
turbine which frequency-modulate, multiplex and transmit sensor outputs; (3) on-site instrumentation to 
record, process and display the sensor output; and (4) statistical analysis of data at the NASA-Lewis 
Research Center in Cleveland, Ohio. 

Two examples of the capabilities of these systems are presented. The first illustrates the 
standardized format for application of statistical analysis to each directly measured parameter. The 
second shows the use of a model to estimate the variability of the rotor thrust loading, which is a 
derived parameter. 


INTRODUCTION 

The U.S. Government has established a Wind 

"Energy Program within the Department of Energy 
(DOE) to encourage the development and promote 
commercial izat ion of wind energy systems. One 
-phase of this program is being managed by the 
"NASA Lewis Research Center (LeRC). An agreement 
“wTth “DOE stipulates that LeRC shall manage both 
the Technology and Engineering Development for 
large ( 100 kW) horizontal axis wind turbines 

(ref. 1). Four wind turbine projects, 

designated the Mod-0 (ref. 2), Mod-OA (ref. 3), 

Mod-1 (ref. 4), and Mod-2 (ref. 5), are part of 
the current development program. In addition to 
these projects, efforts aimed at achieving lower 
machine costs have been initiated. These 
include advanced multi-megawatt (Mod-5) and 200 
to 500 kW (Mod-6) wind turbine projects. With 
regard to all these machines, LeRC has the 
responslbil tty to monitor, and report on their 

performance. 

Despite the diverse characteristics of these 
machines, we can nonetheless identify four 
distinct information/user categories that are 
common to all wind turbines; namely, 
operat ions/the wind turbine itself; 
tecnnology/f ield operations personnel; 

engineering/system and component designers; 
performance/dispatchers or Investors. The 
requirements for each category, as seen in 
FTgure 1, are sufficiently disjoint that 

separate data systems have evolved to meet each 
need. 

At one extreme, with the highest sampling rates, 
are the computer based control systems which 
govern the routine operation of each wind 
turbine generator. The data portion of the 
control system provides information regarding 
adequacy of the wind, status of all critical 


systems, machine alignment with the wind, etc., 
and often monitors more than 100 sensors. For 
certain critical parameters the control system 
must be capable of responding within 
milliseconds. This system is considered to be 
an integral part of the wind turbine and varies 
significantly from one design to the next. 
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Figure 1 - The approximate relationship between 
the number of sensors, characteristic 
sampling time and the function of a 
data system. 
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At the other extreme, with the lowest sampling 
rate, is the performance system to provide 
information for evaluation of a wind turbine's 
availability, reliability and energy 
production. These data requirements are 
generally limited to meteorological, electrical 
and structural parameters with time scales 
ranging from one hour to the lifetime of the 
machine. 

The two remaining information systems. 
Technology and Engineering, are discussed in 
greater detail in the remainder of this report. 
The next section deals with signal conditioning, 
data acquisition and display and the subsequent 
routine digital pre-processing. The third 
section provides an example of the standardized 
statistical analyses. The fourth section 
presents a technique to estimate variability in 
the rotor thrust, which itself is not measured 
directly but is derived from the flatwise 
bending moment. 

TECNOLOGY AND ENGINEERING INFORMATION SYSTEMS 

The Technology Information System has three 
functionally (and spatially) distinct 
components, as seen in Figure 2(a), The data 
signals from the sensors move through the system 
in the sequence: signal conditioning, 
acquisi tion/display, and statistical analysis. 
Physically, these three functions occur: on the 
wind turbine, at or near the base of the wind 
turbine tower and at LeRC, respectively. As was 
seen in Figure 1(a), the Engineering Information 
System is a subset of the Technology Information 
System in that fewer sensors are monitored. In 
addition, the Engineering System uses 
significantly less electronic equipment at each 
wind turbine site. As a consequence, if a 
statistical analysis is desired, some further 
processing and display is performed at the Plum 
Brook Station as illustrated in Figure 2(b). 

The wind turbines monitored by LeRC vary 
considerably in their size, location, blade 
composition and design. Despite, and to some 
extent because of, this variability. It was 
decided that all data system implementations 
should have the same (or functionally 
equivalent) hardware and software. 

Signal Conditioning 

Signal conditioning is performed by a Remote 
Multiplexing Unit (RMU) (ref. 6). As input, an 
RMU can accept up to 32 low-level or high-level 
data signals from a variety of transducers. 
Each RMU contains reference junctions for 
thermocouples as well as the necessary 
electronics for excitation and bridge completion 
of strain gauges. Each DOE/NASA wind turbine 
has one RMU located in the hub, another in the 
nacelle and a third unit at the base of the 
tower in the control room. A listing of all the 
transducers monitored during the initial testing 
of Mod-1 has been reported elsewhere. (ref . 7) 

After a signal is received at the RMU it is 
conditioned (scale and/or offset, amplification 
or attenuation) to a common range and frequency 
multiplexed for output. Each RMU can generate 


two multiplex groups. Each multiplex group 
consists of up to 16 FM signals (+ 125 Hz 
centered at 500 Hz intervals from lOOU Hz thru 
8500 Hz) plus a precise reference tone at 9500 
Hz. Other significant features of the RMU 
include a 4-pole active, Butterworth low-pass 
filter, end-to-end system calibration capability 
(upon command from an external source) and one 
kilometer signal transmission via coaxial cable. 



2(a) Technology Information System 



' FM MULTIPLEX 
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2(b) Engineering inrormation System 


Figure i - Schematic representation of the 

DOE/NASA Wind Energy Technology and 
Engineering Information Systems. 
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Acquisition/Display 

LeRt has two functionally equivalent Technology 
Acquisltion/Display units. One unit is 
installetj in a mobile van and is used to support 
field engineers, at wind turbine sites, thru 
assembly, check-out and initial operation of the 
first unit of each new wind turbine design. It 
is presently located at the Mod-2 site in 
Gcodnoe Hills, Washington. The other unit is 
installed at the LeRC Plum Brook Station, 
Sandusky, Ohio. In addition to processing analog 
tapes received from those wind turbine sites 
equipped with Engineering Acquisition/Display 
units, it is also used in the Supporting 
Research and Technology program based on the 
Mod-0 machine. 

A schematic representation of the electronic 
data processing capability of the mobile unit is 
shewn in Figure 3{a). All RMU-generated FM 
multiplexes entering this unit are, with the 
addition of a time code, recorded in direct 
analog form. This recording capability is 
independent of any other equipment or processing 
activity within the unit. Simultaneously, the 
data can also be routed thru a bank of six sets 
of 16 discriminators which de-multiplex the 
signals and generate analog ( +5V) signals. Any 
all of these 96 analog signals can be 
digitized and routed thru the mini-computer. 
From there they can be processed for real-time 
digital display on the CRT and/or for 

transmittal on digital tape to LeRC for further 
analysis. In addition, any 24 of the 96 analog 
signals may be selected for display on strip 
charts and any single analog signal may be 
routed to a spectrum analyzer for frequency 
content evaluation. All the components shown 
with a stippling in Figure 3(a) can be set up 
and run under computer control at the discretion 
of the unit operator via the control terminal. 
This unit was designed with sufficient capacity 
to simultaneously support up to three wind 
turbines at a single site. 


Even after the initial check-out of a new wind 
turbine design, LeRC retains the responsibility 
to monitor and report the wind turbine’s 
performance. However, the Technology 

Acquisition/Display unit is too elaborate and 
expensive for long term monitoring at each site. 
Therefore, we have identified a functional 
subset, the Engineering Acquisition/Display 
unit, which is installed in the control area of 
each wind turbine. This latter unit provides 
on-line analog display on strip chart (eight 
channels per wind turbine) and continuously 
records 48 signals (as three FM multiplexes) 
plus time code on four tracks of analog magnetic 
tape. The analog magnetic tape recorder 
operates in either of two modes, depending on 
local conditions and requirements. At some 
sites it records until the tape is full (32 
ho\jrs| ' and then automatically rewinds (10 
minutes) and restarts, erasing old data as new 
data are recorded. At other sites the recorder 
operates for 96 hours (by making three passes 
thru the tape using a total of 12 tracks) and 
then automatically turns off. It remains off 
until the tape is replaced and the unit Is 
manually restarted. 


Large volumes of data are often of little value 
in their raw form. Even after processing, they 
may well be of negligible value if the end 
product is overwhelmingly voluminous, 
inadequately disseminated, or excessively 
delayed. To preclude these occurrences we 
routinely perform statistical analyses of both 
the technology and the engineering data. 
Condensed summaries are provided in both graphic 
and tabluar form, using microfiche as the 
distribution medium. 
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3(b) Engineering Acquisition/Display Unit 


Figure 3 - Schematic representation of the 

electronic data processing in the 
Teccnology and Engineering Acquisition/ 
Display Units. In the Technology Unit, 
each device shown with stippling can be 
controlled by the operator from the 
console. 
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Digital , Pre-Processing 

Digital magnetic tapes, can be generated at 
either of the Technology Acquisition/ Display 
units. The digital data consist of 11 readings 
per parameter per (nominal) revolution and are 
stored as a tightly packed, randomly sequenced 
record. Each datum is accompanied by an 
identifying tag. Time markers (to nearest 
millisecond) are merged with the data. 

As the first step on the LeRC main frame system, 
these data tapes are transferred to disk for 
short term storage. During the transfer 
process, the internal representation is 
transformed (in software) from ASCII to EBCDIC. 
The next step compacts this initially large 
dataset ( 5x10^ data values + 5x10^ tags + 

lxl06 time markers) into a more manageable 
form, as follows. The rotor shaft position is 
used to mark the start of each rotation. These 
markers are then combined with the associated 
time markers and processed to give rotor speed 
(rpm) as a function of time. Then, the data 
from approximately 30 sensors of general 
interest are screened to yield maximum and 
minimum values for each parameter for each 
revolution of the rotor. This smaller data set 
( ExlO^data values) is stored on disk and is 
the data base for all further processing. In 
the final step, this latter data set is 
processed onto a microfiche containing the time 
history and statistical summary of each 

parameter of general interest. This entire 

process is shown schematically in Figure 4. 



Figure 4 - Schematic represenatation of the 

final step in information processing 
using the LeRC main computer facility. 


EXAMPLE 1: STANDARD STATISTICAL ANALYSIS 

While the specific set of sensors and their 
associated scale factors will vary from machine 
to machine, nevertheless, a single analysis 
procedure and a common presentation format is 
applied to all data from all machines. 

The data, which are stored in the data base as 
maximum and minimum values for each revolution, 
are transformed to represent the midpoint and 
cyclic values for each revolution. The 
transformation equations are: 

midpoint = (maximum + minimum) / 2 
and 

cyclic = (maximum - minimum) / 2) x 
(1 + f(rpm)) 

where, 

f(rpm) = 2x10-5 (rpm)2. 

The correction function, f(rpm), is introduced 
to compensate for the consistent underestimation 
of cyclic values resulting from the data 
sampling rate of 11 per (nominal) revolution. 

The results of the analysis of each sensor are 
displayed as two frames, one graphical and one 
tabular on the microfiche card. An example of 
the graphical output is shown in Figure 5 and of 
the tabular output in Table I. These data were 
recorded by the on-site Engineering 
Acquisition/Display unit on January 5, 1981, and 
cover more than 12,000 consecutive revolutions 
of continuous machine operation. 

The first frame for each sensor contains three 
graphs. On each graph there are two plots, one 
of the midpoint values with circles as symbols 
and the other of the cyclic values with diamonds 
as symbols. The three graphs are: 

1. Time history. The top graph smooths and 
summarizes the information that one 
associates with a continuous trace on a 
strip chart recorder. One plot is of the 
average, over 30-second intervals, of the 
midpoint values. The other plot is of the 
corresponding cyclic values. 

2. Partitioned distributions (ref. 8). The 
abscissa for this center graph is the wind 
speed midpoint value as measured at hub 
height. These wind speed values are used as 
the basis for sorting the data from the 
sensor of interest. These latter values are 
grouped into subsets such that for each 
subset all the sensor data values were 
obtained at approximately the same measured 
wind speed. Then the data values within 
each subset are separately ranked in 
ascending order. The 16th and 84th 
percentile for each such sequenced subset 
are displayed as horizontal tabs at the end 
of a vertical bar. We also estimate the 
confidence Interval (at the 0.95 level) for 
significant differences of the median^ and 
display it as an interval (denoted by a pair 
of circles or diamonds) on the same vertical 
bar. This entire process is performed 
separately for the midpoint (circle) and 
cyclic (diamond) data values. 
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Figure 5. - Sample frame showing the three graphs produced by 
the standardized statistical processing for each parameter- 
The above graphs comprise over 12,000 rotor revolutions and 
were recorded by the on-site Engineering Acquisition/Display 
Facility. (nt-m = 1.356 ft-lb) . 

Table I. Distributional Information from Figure 5, in Tabular Form 
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3. Cumulative distribution. The bottom graph 
corresponds to a normal distribution, i.e,, 
the abscissa is in units of normalized 
standard deviations and segmented with tick 
marks labelled by percentiles. Such a graph 
has the attribute that if the plotted data 
have a normal (i.e., gaussian) distribution, 
the plot will be linear. For this graph the 
entire set of all midpoint values is 
sequenced and plotted by percentile. This 
process is repeated, separately, for the 
cyclic values. 

The second frame (see Table I) for each sensor 
presents tabular listings of the plotted data 
points from both distributional graphs. Some 
additional (non-plotted) data are also 
tabulated. Because some of these are extreme 
values (i.e., maxima and minima), these latter 
tabulated values must be addressed with caution 
as they might represent spurious noise. 

EXAMPLE 2: ROTOR THRUST STATISTICS 

As an example of using these data systems to 
obtain information about non-measured 
parameters, consider a calculation procedure 
which can be applied to estimate the variability 
in the rotor thrust from the measured flatwise 
bending moment data of a wind turbine blade. 
First an anlysis is made of the changes in the 
flatwise bending moment, which is a measured 
parameter. Then the changes in flatwise bending 
moment are related to changes in rotor thrust, 
which is a derived parameter, and developed as a 
normalized "range" (i.e., variability) 
coefficient. 

Flatwise (Midpoint) Moment Range 

The first step is to produce a pseudo-spectrum 
for the flatwise bending moment. The procedure 
starts with a data set consisting of the 
midpoint values (as defined in the previous 
section and as shown in figure 5) of the 
measured bending moment for each of some large 
number, n, of rotor revolutions. From these 
data, the range, i.e., the maximum minus the 
minimum, is calculated repeatedly over 
successively doubled time intervals. That is to 
say, the calculation of range values is made 
over every two cycles, again over every four 
cycles, again over every eight cycles, etc. 
This procedure yields n range values of the 
midpoint flatwise bending moment. 

These range values are sorted in ascending order 
and plotted as a standard normal probability 
graph, as shown in Figure 6, using the scale 
shown to the left. The data for this particular 
calculation are the same data as seen earlier in 
Figure 5. The flatwise bending moment was 
measured near the root end of the blade at 
approximately five percent of the span, i.e., 
100 cm from the axis of rotation. The median of 
the range of the flatwise bending moment data is 
8,800 N-m and the 84th percentile is 24,400 N-m. 

These same data are replotted in Figure 6 for 
the same abscissa, but using as the ordinate the 
logarithmic scale to the right of the graph. 
The fit of the straight line to these 


data from the 50th to the 99th percentiles 
indicates a log-normal distribution that is 
truncated at the very high end. When the same 
analysis is made on data which include 
start-up/shut-down sequences, there is no 
truncation of the distribution and log-normality 
prevails to the uppermost extent of the data set 
(ref. 9). In either case the results in the 
region from the 50th to 99th percentile are 
essentially the same. Therefore, the 
information extracted is drawn from the central 
region of the distribution which is, in 
statistical jargon, more robust. 

Rotor Thrust 

The conversion from mean flatwise bending moment 
to mean rotor thrust is made by using the 
MOSTAB-WT (ref. 10) simulation program to 
determine the ratio, C, of the change in thrust 
to the change in flatwise bending (N/N-m). For 
wind speeds below the rated wind speed (11 m/s) 
C = 0.016, and for wind speeds above the rated 
wind speed, C = 0.011. Since both conditions 
occurred during the machine operation and the 
calculations described above do not distinguish 
between these two cases, an average of the two 
cases, C = 0.013, is acceptable for this 

analysis. Applying this average value to the 
flatwise moment range data in Figure 6, we 
obtain the average results shown in Table II. 
The median value of the range of the rotor 
thrust is 1600 N and the 84th percentile is 4400 
N, This information can be normalized by the 
rotor thrust at rated wind speed (32,000 N) and 
restated as a dimensionless parameter. This 
gives the dimensionless values of 0.05 at the 
50th percentile and 0.14 at the 84th percentile 
for the ratio of the range of rotor thrust to 
the rotor thrust at rated wind speed. 
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Figure 6 - Normal probability plots for the range 
of variability in the flatwise bending 
moment, based on the same data as in 
Figure 5. 
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Table II Rotor Thrust Varability for Mod-OA 

change in rotor thrust , 

change in flatwise bending moment = 0.013m' 

rotor thrust at rated wind speed= 32,000 nt 


Parameter 

Median 

flatwise moment range 


{midpoint), nt-m 

8,800 

^ctor^ thrust range. 


nf-m 

1600 


rotor thrust range 

rotor thrust at rated wind speed 0.05 


84th Percentile 


24,400 

4400 

0.014 


CONCLUDING REMARKS 

Two related data systems, one for technology 
.development and the other for engineering 
evaluation, have been identified and described. 
Typical data, in this instance six consecutive 
hours of flatwise bending moment data were 
presented. The description of the extraction of 
Information thru the application of statistical 
analyses Included both the standardized analysis 
procedure and an analysis specific to evaluating 
rotor .thrust,. 
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QUESTIONS AND ANSWERS 
H.E. Neustadter 


From: 

Q; 

A: 

From: 

Q: 

A: 

From: 

Q: 

A: 


G . McNerney 

How extensive is the data base for each site, (sample rates, total time, and chan- 
nels sampled) , and is the data available in raw form for dissemination outside NASA? 

The raw form of the data is analog recording of FM encoded data and is not availa- 
t>le for general dis tr'ihu.t'Lon • For d'vg^tal analysts the sampling rate ts 21 ttmes 
per revolution for up to 48 channels over a six hour duration for each analysis 
performed. There are about 100 analyses currently on microfiche and the production 
rate is, nominally, one per machine per week. 

W.N. Sullivan 

What is the turnaround time on creation of engineering information data sets? 

Under ideal circumstances it could he 2 or Z days. Typical turnaround is 7 to 12 
days . 


P.C. Klimas 

1) Does your data system measure rotor C^? If so, is the data available? 

2) What is the time interval used for the measurements which you "bin"? 

1) Not yet, because most sites do not have free-stream wind data available. 

2) The values which are '^binned'^ are the midpoint and cyclic values for a revo- 
lution. The time interval depends on the rotor speed. 
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